Reply to Comment on "Spherical 2+p spin-glass model: an analytically solvable 

model with a glass-to-glass transition" 
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In his Comment, Krakoviack [Phys. Rev. B (2007)] finds that the phase behavior of the s + p 
spin-glass model is different from what proposed by Crisanti and Leuzzi [Phys. Rev. B 73, 014412 
(2006)] if s and p are larger than two and are separated well enough. He proposes a trial picture, 
based on a one step replica symmetry breaking solution, displaying a mode-coupling-like glass-to- 
glass transition line ending in a A3 singularity. However, actually, the physics of these systems 
changes when p — s is large, the instability of which the one step replica symmetry breaking glassy 
phase suffers turns out to be so wide ranging that the whole scenario proposed by Krakoviack must 
be seriously reconsidered. 

PACS numbers: 75.10.Nr, ll.30.Pb, 05.50.+q 
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The model under consideration consists of N spherical 
spins, i.e., continuous variables <Xj ranging from —00 to 
00, obeying the spherical constraint 



iV 



(i) 



and interacting via two different random quenched multi- 
body interactions. The Hamiltonian of the spherical s+p 
spin glass model is, then, defined by 

ii<...<i s ii<...<i p 
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where jf*] ■ ,t=s,p, are uncorrelated zero mean ran- 
dom Gaussian variables of variance 



(j [t) V 

Y i1i2-.it J 



2N*~ ] 



H < 



< it 



(3) 



The scaling with N guarantees a correct thermodynamic 
limit. As one can see from Eq. ([2]), we are considering the 
mean-field approximation, in which each spin interacts 
with all other spins. Notice that only distinct s-uples 
and p-uples are taken into account in the Hamiltonian. 

The properties of the spherical s+p model strongly 
depend on the values of s and p: for s — 2, p = 3 the 
model reduces to the usual spherical p-spin model in a 
fields with a low temperature one step Replica Symme- 
try Breaking (RSB) phase (i.e., the "mean-field glass" 
phase), while for s — 2, p > 4 the model possesses an 
additional Full RSB low-temperature phased The model 
case under investigation will be here the one with both s 
and p larger than 2. 

One of the interesting features of the model is that by 
defining the auxiliary thermodynamic parameters fi p = 
p0 2 Jp/2 a straightforward connection can be made with 
the mode-coupling theory (MCT). In the high tempera- 
ture regime (i.e., in the fluid phase), indeed, the dynamic 
equations of the model can be formally rewritten as MCT 



equations in terms of the /i's and an exact mapping can 
be set with mode coupling schematic theories F s —i p—i 
with a scalar kernel. In particular, the F13 theory stud- 
ied by Gotze and Sjogren^ is dynamically equivalent to 
a 2 + 4 spherical spin model. 

A partial analysis of the phase space of the 2 + 4 model 
was carried out in Ref. [H where, however, only the dy- 
namical stability of the 1RSB phase was considered leav- 
ing out a large part of the phase space and, in particular, 
the question of the transition between the 1RSB and the 
FRSB phases. In that analysis, indications for a glass- 
to-glass (G-G) transition line, called A2 line in MCT, 
ending in a A3 point (namely, a cusp) was found. How- 
ever, the whole line was shown to remain in a region of 
the parameter space where the 1RSB phase was proven 
unstable. The study of the dynamics and of the statics of 
the generic spherical 2+p spin model has been completed 
in further works by the authors^i^, especially showing 
that a different "glassy" phase arises (the "1-Full" RSB 
phase), different from the 1RSB one, as well as a Full 
RSB phase, typica l of systems with discrete spins (see, 
e.g., Refs. I9l lldllllfl2l ). In Fig. [T]we reproduce the phase 
diagram of the 2 + 4 model 7 to illustrate this. 

If both s and p are larger than 2, it was shown in Ref. 
7 that no 1-FRSB, nor FRSB phases occur and one has 
a smooth transition between a fluid/paramagnetic (PM) 
and a glassy phase (1RSB) (first dynamic, then static, 
lowering the temperature). Krakoviack 1 pointed out in 
his Comment that this is true only as far as p — s is 
not too large. Indeed, he analyzes a parameter region 
of the spherical s + p-spin model not considered before, 
where the difference between s and p is large (in particu- 
lar s = 3, p = 16) providing hints for a transition between 
two glassy phases both presumably described by a 1RSB 
solution. A transition line of the type A2, was drawn 
in the phase diagram terminating in a A3 points In the 
RSB theory language, this corresponds to a point where 
a cusp between a physical and a non-physical branch of 
the constant m = 1-line occurs, where m is the value of 
the RSB parameter x at which the overlap step function 
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FIG. 1: Phase diagram of the spherical 2 + 4 spin glass 
model. RS: replica symmetric (paramagnetic) phase; 1RSB: 
one-replica symmetry breaking phase; FRSB: full replica sym- 
metry breaking phase; 1-FRSB: one-full replica symmetry 
breaking. The dashed lines refer to dynamic transitions. The 
continuous transition between the RS and the FRSB phases 
and between the FRSB and 1-FRSB phases are the same for 
statics and dynamics. Inset: close-up of the region around 
the "end point" . 
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FIG. 2: Plot of the constant m line at m = 1, cf. Eqs. 
0-©, in the /lie — /13 phase diagram. This is the dynamic 
transition line between the paramagnetic (PM) and the glassy 
1RSB phases. The portion between the points T and A3 is 
called A2 line in MCT. Following the arrows from /116 = the 
complexity first decreases below zero and then increases. The 
dashed part of the curves correspond to points of the phase 
diagram for which the complexity is negative. 



q{x) discontimiously changes value. This behavior is very 
reminiscent of what happens in the above-mentioned 2+p 
spherical spin model. 

The dynamic transition line (constant to = 1) is plot- 
ted in Fig. [2 to be compared with Fig. 1 of Ref. QJ It is 
determined by solving the equations 
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(a - l)gi - (s - 2) 
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For large enough p — s the line displays, e.g., for s = 
3,p = 16 (Fig. [2]), the so-called swallowtail. 

Let us now consider in more detail this dynamic tran- 
sition line. In particular, we concentrate on the vertical 
part of the line above the crossing point: the A2 line TA3, 
cf. Fig. [2] This is the first of the G-G lines proposed 
by Krakoviack (Fig. 1 of Ref. [J). In Fig. [Hit is shown 
that, above a given fj, s value, the solutions along the line 
correspond to states of negative complexity, excluding, 
therefore, the A3 point and reducing the length of the 
candidate G-G transition line (to an upper bound that 
we denote by A3, where the complexity is zero). 

Analyzing the model further on, one can actually ob- 
serve that, in the model cases where the swallowtail and 
A2 G-G line show up (e.g., for s = 3 it occurs when 
P > 10), in some region of the phase diagram the hy- 
pothesized 1RSB solution turns out to be inconsistent 
and must be, therefore, substituted by different, more 
complicated Ansdtze. One sufficient signature of this in- 
stability, by no means necessary, is, e.g., that the com- 
plexity of the system becomes negative. We show in Fig. 



H]the constant m = 1-line, i.e., the dynamic PM/1RSB 
transition line. Starting from [Xiq = (and following the 
arrows) the complexity of the system first decreases from 
positive to negative and, afterwords, increases along the 
vertical branch, becoming positive again. Part of the A2 
line thus corresponds to a system with negative complex- 
ity. 

A serious danger exists that the instability might ex- 
tend into the frozen region, heavily affecting the G-G line, 
similarly to what occurred in the case s = Deepen- 
ing the analysis of the model one finds out that this is 
actually what happens. 

To show how the instability occurs, we first begin with 
the 1RSB solution, looking at its "bugs". In Fig. [3] the 
loci of zero complexity are reported. These lines are con- 
structed following the m-lines with to < 1 as we did above 
for the 1-line. For small to the m-lines are single-valued 
in the (file, fj.3) plane and the complexity S computed 
along an m-line first becomes negative and then positive 
as //i6 increases. Increasing to, the TO-lines start to be 
multivalued in some parameter region (see the curve for 
to = 0.6 in Fig. [3]). Moreover, as to further increases to 
to < 1, the generic m-line crosses itself and, always start- 
ing from /ii6 = 0, the point at which £ first decays to zero 
turns out to be on the right side with respect to the point 
where S reaches zero from below (look, e.g., at the inset 
of Fig. [3|). Spanning the plane with m-lines, two lines 
of loci of zero complexity system points are constructed, 
signaling that the 1RSB solution is not the physical one 
and must be rejected and substituted. 

Apart from the dynamic A2 line examined above, 
Krakoviack 1 puts forward also other putative G-G tran- 
sition lines, both in statics and in dynamics. Without 
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FIG. 3: Zero complexity curves and constant m-lines. Three 
m-lines are explicitly plotted as examples: m = 1 (dynamic 
transition), m = 0.6 and m = 0.84 (in the inset). As a guide 
for the eye we have put arrows: starting from /ii6 = 0, E 
decreases from positive to negative values and, then, increases 
again to positive values. With respect to the E = curve, 
the side of the phase diagram where the points of positive 
and negative complexity are, depends on the value of m. As 
m increases towards one, the m-lines become more and more 
entangled. The full curves represent the parts of the m-lines 
whose points have E > 0, the dashed curves represent the 
points with E < 0. The two zero-complexity curves (dashed- 
dotted lines) start from the 1-line (empty diamonds). 




10 20 30 40 50 60 70 80 

he 

FIG. 4: Static phase diagram of the 3 + 16 model in the (fiie, 
fis) plane in the parameter region considered by Krakoviack 
in Fig. 2 of Ref. 1. His putative G-G transition line (TC) is 
also plotted: the whole line is contained in the 2RSB thermo- 
dynamically stable glassy phase. 



entering in a critical analysis of their derivation, we re- 
call that the static line (Fig. 2 of Ref. [l]) corresponds 
to the line where the free energy of two different (but 
both 1RSB) glassy phases is equal, that is to a first order 
phase transition. Looking for RSB schemes of computa- 
tion yielding a stable thermodynamics everywhere in the 
parameter space for large p — s, one observes, however, 
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FIG. 5: Dynamic phase diagram of the 3 + 16 model in the 
(M3, Mie) plane to be compared with Fig. 3 o^. The AB line 
corresponds to the A2 transition line terminating in B=A° 
that we consider in Fig. [2] and is not a true transition line. 
The AC line is the dynamic analogue of the TC line in the 
static case, cf. Fig. [4] 



that a qualitatively different phase appears in the region 
of the phase diagram where this G-G transition line runs: 
a two step RSB stable phase. We reproduce Krakoviack's 
static G-G line in Fig. [U showing how it is "eaten" by 
the 2RSB phase. 

The complete stabilized scenario for the spherical s+p 
model (with s,p > 2) is, actually, even more complicated 
than in the 2+p case and leads to new results in spin-glass 
theory and provides the mean-field analogues of glasses 
with Johari-Goldstein processes^ (i.e., thermalized pro- 
cesses with large relaxation times, though shorter than 
the structural relaxation time), cf. Ref. Il4l Here we limit 
ourselves to show, in Fig. [4] the phase diagram computed 
in the framework of RSB theory around Krakoviack's pu- 
tative static G-G line and to remark that the TC line is 
ruled out by the presence of the 2RSB phase. 

In Fig. [5]we show the diagram for the dynamic tran- 
sition lines. The AC line is the dynamic equivalent of 
the static TC line of Fig. [5J that is, the loci where 
two apart (1RSB) solutions display the same complexity, 
rather than the same free energy. We use for the end- 
points the same notation used in Ref. Q], Fig. 3, where B 
is what we called A3 point above. The candidate transi- 
tion line AC turns out to be embedded in the 2RSB phase 
(the 1RSB/-2RSB and 2RSB-1RSB// transitions are here 
dynamic transitions) and is, thus, discarded. The candi- 
date AB line, instead, is not even a transition line, since 
the complexity of each of its points is lower than the com- 
plexity along any of the crossing (horizontal) m-lines. 

Although the putative G-G lines proposed by Krakovi- 
ack for the present model and, in particular, the anal- 
ogy with the G-G transitions conjectured in the frame- 
work of MCT are ruled out^ we stress that this does 
not mean that transitions between qualitatively different 
glassy phases are absent in the spherical s+p spin model 
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with s and p larger than two and large p — s. See, e.g., 
Figs. [41151 

In summary, Krakoviack's observation that if p — s is 
large enough G-G transitions show up is valid. However, 
comparing Figs. 1, 2 and 3 of his Comment with, respec- 



tively, Figs. [2l |4] and [5] in this paper, one can conclude 
that those G-G transitions proposed in the Comment, 
and the consequent supposed scenario, have to be re- 
jected, even from a heuristic point of view. 
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